Toward TeV Conformality 
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We study the chiral properties of an SU(3) gauge theory with Nf massless Dirac fermions in the fun- 
damental representation when Nf is increased from 2 to 6. For Nt = 2, our lattice simulations lead to a 
value of (V>V>) / F 3 > where F is the Nambu-Goldstone-boson decay constant and (ipip) is the chiral conden- 
sate, which agrees with the measured QCD value. For Nf — 6, this ratio shows significant enhancement, 
presaging an even larger enhancement anticipated as Nf increases further, toward the critical value for 
transition from confinement to infrared conformality. 



PACS numbers: ll.IO.Hi, ll.15.Ha, ll.25.Hf, 12.6Q.Nz 

Introduction Theories with an approximate conformal 
symmetry could play a role in describing new physics at 
the TeV scale and beyond. While a nonsupersymmetric, 
vectorlike gauge theory exhibits confinement and sponta- 
neous chiral symmetry breaking with a small number Nf 
of massless fermions, it becomes conformal in the infrared, 
governed by a weak infrared fixed point if Nf is larger, but 
just below the value for which asymptotic freedom sets in 
HI]. There is evidence from lattice simulations 0-01 that 
this infrared conformality persists down through a "con- 
formal window" of Nf values where the fixed point can 
become strong, and that a transition to the confining and 
chirally broken phase takes place at some value Nf. 

Even for Nf < Nf there can remain an approximate in- 
frared fixed point providing that < Nf — Nf <C A^. 
The scale F of chiral symmetry breaking is then small rel- 
ative to the intrinsic scales of the theory, and the fixed point 
approximately governs the theory from the breaking scale 
out to some higher scale. 

This "walking" phenomenon can play an important 
phenomenological role in a technicolor theory of elec- 
troweak symmetry breaking. Flavor-changing neutral cur- 
rents (FCNC's), which are present when the technicolor 
theory is extended to provide for the generation of quark 
masses, can be too large unless the associated scale K ET q 
is high enough. But then the first- and second-generation 
quark masses are typically much too small. They are pro- 
portional to the quantity (ipxp) /A 2 ETC , where tp is a tech- 
nifermion field and {ipip) is the bilinear fermion conden- 
sate defined (cut off) at A ETC - Walking can lift the quark 
masses by enhancing the condensate (ipi/j) significantly 
above its value [0(47rF 3 )] in a QCD-like theory 
while keeping A% TC large enough to suppress FCNC's. 



The enhancement of (tpip) / F 3 as Nf — > Nf from be- 
low has been indicated by Feynman-graph-based studies. 
But it is important also to use nonperturbative methods 
since the couplings involved are strong. This Letter de- 
scribes a first step in this program. We focus on an SU (3) 
gauge theory with Nf massless Dirac fermions in the fun- 
damental representation. Lattice studies have shown that 
the Nf = 8 theory is chirally broken, with no evidence 
for even an approximate infrared fixed point J2, EL 01; 
while there is lattice evidence for conformal behavior at 
N f = 12, indicating that 8 < N c f < 12 0-@]. 



We present results here for the values N 



/ 



2 and 



Nf = 6, drawing on newly available computational re- 
sources, including 150 x 10 6 core hours on the BlueGene/L 
supercomputer at Lawrence Livermore National Labora- 
tory (LLNL). Starting with Nf = 2 allows us to check 
the reliability of our methods by comparison with the phe- 
nomenological value of (^ip)/F 3 for QCD. Proceeding 
carefully toward Nf is prudent since the emergence of 
widely separated scales associated with the approximate 
infrared fixed point of walking is problematic for lattice 
methods. 

Methods For a range of small fermion masses m, we 
compute the Nambu-Goldstone-boson (NGB) mass M m , 
the NGB decay constant F m , and the chiral condensate per 
fermion (ipij)) m . To set a physical scale, we also com- 
pute the mass M p m of the analogue of the p meson and 
the Sommer scale r 0jTO at which r 2 dV(r)/dr = 1.65, 
where V(r) is the static potential II 1 ill - Since our goal is 
to search for the enhancement of {ipijj) /F 3 as Nf —> Nf, 
from the emergence of walking between the physical length 
scale and the ultraviolet cutoff, taken here to be the lat- 
tice spacing, it is important to keep the lattice spacing fixed 



2 



(and small) in physical units. We first choose a value for 
(3 = 6/g 2 at Nf = 6, giving a physical scale of several 
lattice units. For Nt = 2, we then tune (3 to match the 
same physical scale in lattice units. 

For small enough fermion mass (and yet large enough 
to insure that finite- volume effects are small), the extrap- 
olation m — > can be carried out by fitting the results 
for M^, F rn and (ipip) m to continuum yPT. The next-to- 
leading-order (NLO) expressions are fl2ll 
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where z = 2{ip^p) / (Ait) 2 F 4 . The leading terms in- 
corporate the Gell-Mann-Oakes-Renner (GMOR) relation. 
These expressions will be directly useful for Nf = 2, but 
because of the growth with Nf of the chiral-log terms, not 
for Nf = 6. 

Simulation Details We use domain-wall fermions with 
the Iwasaki improved gauge action, as used by the RBC- 
UKQCD collaboration lfl3ll . Lattice fermion discretization 
typically breaks chiral symmetry, but in the domain-wall 
formulation the breaking is exponentially suppressed (with 
flavor symmetry preserved), making it ideal for the study 
of chiral dynamics. Gauge configurations are generated us- 
ing the hybrid Monte Carlo method as implemented in the 
USQCD application libraries, in particular CPS, via a mul- 
tilevel symplectic integrator, and using Hasenbusch pre- 
conditioning and chronological inversion. Autocorrelation 
is reduced by blocking over sets of 50 trajectories. 

The lattice volume is set to 32 3 x 64, with the length of 
the fifth dimension L s = 16 and the domain-wall height 
m = 1.8. All quantities are given in lattice units. For 
Nf = 6 we choose j3 = 2.10. For Nf = 2 the choice 
f3 = 2.70 then leads to nearly the same physical scale in 
lattice units. Simulations are performed for fermion masses 
wij = 0.005 to 0.03. At finite lattice spacing, even with 
wij = 0, the chiral symmetry is not exact, with the vi- 
olation captured in a residual mass m res <C mf. For 
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2.60 x 10~ 5 , while for N 
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The total fermion mass m is then 



^, TTl res 

m res = 8.23 x 
m = mf + m res . 

Although global topological charge, Q, is an irrelevant 
quantity with infinite volume, in a finite volume it becomes 
relevant IU4I1 . On a discretized lattice, Q is not conserved, 
with the system evolving between sectors, an evolution cru- 
cial for the correct sampling of the path integral at finite 
volume. With very light fermions, the evolution of Q slows 
dramatically using current Monte Carlo methods 11511 . We 



find that Q evolves sufficiently for mf > 0.01 for both 
Nf = 2 and 6. At mf = 0.005 it does not, leading to sys- 
tematic shifts in (^V')m an d F m , which we will explore in 
a future paper. Here, we present results for mf = 0.005, 
but do not include them in our analysis. This also ensures 
that for each m, M m L > 4, keeping the NGB Compton 
wavelength well inside the lattice (the p regime). Results 
for mf = 0.025, 0.03 are likewise not used in our analysis. 

Results We first report our results for the extrapolated 
values of the physical scales M p m and l/r , m , with /3 = 
2.10 at N f = 6 and (3 = 2.70 at N f = 2. For small 
enough m, they can be extrapolated torn = using yPT 
fl^l . where the NLO terms are now linear in m (there is 
no m log(m) term). This and the small change in these 
quantities in the range m = 0.01 - 0.02, indicate that a 
linear extrapolation should suffice for both Nf = 2 and 
Nf = 6. The extrapolated values in lattice units are r^ 1 = 
0.111(4) (N f = 2) and 0.100(6) (N f = 6), and M p = 
0.198(14) (Nf = 2) and 0.207(15) (Nf = 6). Thus, to 
within the 10% accuracy of this Letter, both r^ 1 and M p 
remain fixed going from Nf = 2 (with j3 = 2.70) to Nf = 
6 (with j3 = 2.10V We note finally that for QCD, r = 
0.378(9) GeV" 1 IH, giving 1/M p r = 0.488(12), in 
reasonable agreement with our Nf = 2 value l/M p r Q = 
0.561(44). 

We now turn to M^, F m , and (iptp} m , noting first 
that M 2 n j2mF m extrapolates to (tjj^/F 3 (the GMOR 
relation) in the chiral limit. We can get an estimate of 
the enhancement of this ratio by comparing M^ l /2mF rn 
for Nf = 6 to that for Nf = 2, at finite m f . 
We do this by plotting the ratio of ratios R m = 
[Ml/2mF m ] 6 f/[Ml/2mF m ] 2f in Fig. ffl We use the 
ratio M 2 n /2mF m here, but we could also use the ratios 
$4>} m /F r l or (M, 2 „/2m) 3 / 2 /(^V) I 1 ,{ 2 : which also ex- 
trapolate to {ijjip)/F z , and which show the same trend. 
The evident trend in Fig. Q] is that R rn increases as m/ 
decreases. Even disregarding the point at nif = 0.005, 
this suggests that the extrapolated value will be well above 
unity. A linear extrapolation of the data for m = 0.01 — 
0.02 gives 1.58(19) for m = 0, implying an enhancement 
in this range or above unless there is a downturn in R m . 
This would require either special values of the x?T param- 
eters given the natural upturn of the combined chiral-log 
terms in Eqs. Q] and 12 or a significant downturn before 
%PT turns the curve up again as m — > 0. 

The separate simulation results for M? n /2m, F m , and 
(ipvjj) m are shown in Figs. [2-E1 We discuss first the ev- 
idence that for Nf = 2, NLO x?T (Eqs. £Q - ®) can 
be used for the extrapolation to m = 0. We note that for 
the range m f = 0.01 - 0.02, for which M m L > 4, we 
have F m L = 0(1), raising another concern about the ap- 
plicability of NLO xPT- Nevertheless, it has been observed 
that for QCD studies with similar F m L values, the finite- 
volume corrections are no more than a few percent IU3I1 . 
We will examine these finite-volume corrections in a fu- 
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FIG. 1: R m = [M^/2mF m ] 6f /[M^ l /2mF m ] 2 f, versus m = 
(ra 2 f + m e j)/2, indicating enhancement of (ipip) / F 3 at Nf — 6 
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ture paper [18], relying here on the smallness of the RBC 
finite- volume corrections . 

From Figs. |2]and[3l we see that for N f = 2, M^/2m 
and F m change little in the range irif = 0.01 — 0.02, in- 
dicating along with the Nf dependence in Eq. [3]that NLO 
%PT should provide a reliable fit. Our results for (ipip) m 
show the expected dominance of the large linear contact 
term (which does not preclude the use of xPT). We there- 
fore carry out the Nf = 2 extrapolation to m = using 
the combined NLO chiral expansions of Eqs. ([]]) - (0. The 
five-parameter fit, shown in Figs. [2] - |4l leads to extrapo- 
lated values F = 0.0209(41) and {^}/F 2 = 0.99(17), 
giving the NLO xPT result 
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{W>) 

F 3 



47.1(17.6). 



(4) 



The other fit parameters are olm = 0.31(62), ap = 
0.64(47) and a c = 83(29). The values of the fit param- 
eters excluding ac indicate that we are near the limit of 
applicability of x?T- We note also that x 2 /d.o.f. = 6.50 
with an uncorrected fit and 4 degrees of freedom. 

We next compare our Nf = 2 results for {^tp)/F 3 
and M p /F to the QCD quantities (qq)/f% and m p / f„, 
a reasonable comparison since the light-quark masses m q 
are so small. With / w = 92.4(0.3) MeV and m p = 
775 MeV, we have m p /f„ = 8.39(0.04), compared 
to our value M p /F = 9.4(2.5). The condensate (qq) 
is renormalization-scheme dependent, as is m q . In the 
MS scheme at 2 GeV (~ 2.6mA Ref. [19] finds 



{QQ)2 Gev/fn = 24.1(4.3). In our case, (tptp) is defined 
by lattice regularization with a -1 ~ 5M P (equivalent to 
3.85 GeV). The increase in {qq) going to this higher scale 
can be estimated perturbatively from the anomalous dimen- 
sion of the mass operator [20]. We find ((7(7)3.85 Gev//^ = 
29.5(5.3). There is also a renormalization factor Z MS , 
which converts the MS condensate to the lattice-cutoff 
scheme. Using Ref. JH, we find Z MS (3.85 GeV) = 



FIG. 2: The slope of the pseudoscalar mass squared M m /2m in 
lattice units, as a function of fermion mass. The fit for Nf — 2 
is a joint fit to M m , F m and (i/)t/i) m , using the (solid) points at 
mf = 0.01 - 0.02, constrained to match NLO xPT. 
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FIG. 3: The Goldstone-boson decay constant F m in lattice units, 
as a function of fermion mass. The fit for Nf = 2 is a joint 
fit to M^, F m and (tp^m, using the (solid) points at mj = 
0.01 - 0.02, constrained to match NLO xPT. 
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FIG. 4: The chiral condensate per fermion (ipip) m in lattice units, 
as a function of fermion mass. The fit for Nf = 2 is a joint 
fit to M m , F m and (tfjift) m , using the (solid) points at mj = 
0.01 - 0.02, constrained to match NLO xPT. 
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1.227(11), and therefore (qq} 3 . 85 GeVM /f 3 = 36.2(6.5), 
which agrees with Eq.|4] within errors. 

Finally, we discuss our Nf = 6 results, where both the 
simulation data and the Nf dependence of the NLO-^PT 
expressions indicate that the m values are not yet small 
enough, and the volume not yet large enough, to validate 
the use of NLO xPT- (We nevertheless note that a sim- 
ple polynomial fit to the solid points for M 2 n /2m, F m , 
and (ipijj) m leads to extrapolated values that satisfy the 
GMOR relation at 1.3a.) We argue, though, that a conser- 
vative lower bound can be indirectly placed on {ipip) /F 3 
for Nf = 6 by bounding F from above and (tpip)/F 2 
from below. 

For F rn (Fig. g), the N f = 6 points at m f = 0.01-0.02 
decline steeply with decreasing m. The eventual reliabil- 
ity of [Eq. [2) at lower masses will, because of the 
negative curvature in the chiral-log term, bend the points 
down even more rapidly. An upper bound on the extrap- 
olated value F should therefore emerge from a linear fit 
through the three points. This gives F < 0.0208(26), es- 
sentially the same as the value of F in the Nf = 2 case. 
For M^/2m (Fig. O, the Nf = 6 points in the range 
nif = 0.01 — 0.02 are nearly flat as a function of m. Since 
XPT behavior [Eq{TJ, with its positive curvature, sets in ei- 
ther in this mass range or lower, bending the points up as 
m is decreased, a lower bound on the extrapolated value 
(tp^)/F 2 should emerge from a linear fit through these 
three points. This gives (vjjip)/F 2 > 1.25(5). 

Together, these give the conservative lower bound 

N f = 6: > 60.0 (8.0). (5) 

The central value of even this conservative bound shows 
roughly a 30% enhancement relative to the Nf = 2 central 
value [Eqj4j. A comparison to the more precise Nf = 2 
ratio from QCD phenomenology [36.2 (6.5)], leads to an 
increase of at least 30% at la. And we note that these 
bounds are well below the enhancement indicated by in- 
spection of Fig. Q] 

Conclusion The ratio (tpip)/F 3 in an SU(3) gauge 
theory with Nf massless Dirac fermions in the fundamen- 
tal representation, with simulations carried out with a small 
lattice cutoff fixed in physical units, is enhanced when Nf 
is increased from 2 to 6 - by 40% or more from inspec- 
tion of the simulation results (Fig. Q]). An enhancement of 
less than 40% would require a significant downturn in R m 
(Fig. [Q before the eventual applicability of NLO yPT very 
likely bends it up again. Even the conservative lower bound 
of Eq. [Syndicates a substantial increase. This enhancement 
of (ipip) I F 3 at Nf = 6 relative to Nf = 2, arising from 
distance scales between the confinement scale (~ 1/M p ) 
and the lattice scale (~ (5M P ) _1 ), appears even though 
the Nf = 6 theory is not yet walking. (The running cou- 
pling at the lattice scale for Nf = 6 is still rather weak - 
not much stronger than for Nf = 2.) 



For comparison, a perturbative MS computation of the 
enhancement of {ipip) / F 3 at Nf = 6 relative to Nf = 2, 
based on the anomalous dimension of {ipip), gives an en- 
hancement on the order of 1.05 to 1.1, depending on the 
order of perturbation theory ifTill . This can be converted to 
the lattice scheme using the results of Ref. 112 ill by multi- 
plying by the ratio Z^f/Zff = 1.449(29)/1. 227(11) = 
1.18(3). The result is a perturbative enhancement on the 
order of 20% - 30%. 

It will be helpful to obtain results for smaller m (with a 
larger lattice volume) and perhaps study the chiral extrap- 
olations at NNLO. We are now exploring larger values of 
Nf (— > Nf) and other gauge groups, and we will study the 
consequences of walking for quark and lepton mass gener- 
ation and electroweak precision measurements. 
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